Load data: sintax = 0.6
Load data: sintax = 0.8
# Functions:
source(paste0(WD, "scripts/finding_the_core.R"))
source(paste0(WD, "../2023_RETHINK_sample_collection/scripts/functions/PCOA_from_long_data.R"))
source(paste0(WD, "../2023_MiDAS_OK_testing/script/ADONIS_from_long_data.R"))
# Make PCoA on the differnet taxonomic levels
source(paste0(WD, "../2023_RETHINK_sample_collection/scripts/functions/PCOA_from_long_data.R"))
source(paste0(WD, "scripts/functions/growing_dying_bar_plot.R"))
SampleContent2_colors <- c("#1e8449","#21618C","#b03a2e", "#ffb74d", "#e65100", "#795548")
# Functional groups
Nitrifiers = c("g__Nitrosomonas","g__Nitrosospira","g__Nitrospira","g__Nitrotoga")
Denitrifiers = c("g__Zoogloea", "g__Rhodoferax", "g__Thauera", "g__Rhodobacter",
"g__Sulfuritalea", "g__Paracoccus", "g__Azoarcus")
PAOs <- c(
"g__Ca_Phosphoribacter",
"g__Azonexus",
"g__Ca_Accumulibacter",
"g__Ca_Lutibacillus",
"g__Tetrasphaera",
"g__Microlunatus" )
GAOs <- c("g__Ca_Competibacter",
"g__Defluviicoccus",
"g__Propionivibrio",
"g__Micropruina",
"g__Ca_Contendobacter",
"g__Ca_Proximibacter")
Filaments <- c("g__Ca_Microthrix", "g__Leptothrix","g__Sphaerotilus","g__Ca_Villigracilis","g__Trichococcus",
"g__Thiothrix","g__Ca_Promineofilum","Haliscomenobacter","g__Gordonia","g__Sarcinithrix",
"g__Ca_Amarolinea","g__Kouleothrix","g__Ca_Alysiosphaera",
"g__Nocardioides","g__midas_g_1668","g__Anaerolinea","g__Ca_Caldilinea",
"g__Ca_Hadersleviella","g__midas_g_344",
"g__Skermania","g__Ca_Nostocoida","g__Neomegalonema","g__Beggiatoa",
"g__Ca_Brevefilum")
df_functional <-
tibble(Genus = c(Nitrifiers, Denitrifiers, PAOs, GAOs, Filaments),
functional_group = c(rep("Nitrifiers", length(Nitrifiers)),
rep("Denitrifiers", length(Denitrifiers)),
rep("PAOs", length(PAOs)),
rep("GAOs", length(GAOs)),
rep("Filaments", length(Filaments))
))
NOT MADE YES